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Abstract 

Neutrinos with hierarchical masses and t,wo large raixing angles may naturally originate 
from sequential dominance (SD). Within this framework we present analytic expressions 
for the neutrino mixing angles including the next-to-leading order (NLO) and next-to-next- 
to-leading order (NNLO) corrections arising from the second lightest and lightest neutrino 
masses. The analytic results for neutrino mixing angles in SD presented here, including the 
NLO and NNLO corrections, arc applicable to a wide class of models and may provide useful 
insights when confronting the models with data from high precision neutrino experiments. 
We also point out that for special cases of SD corresponding to form dominance (FD) the 
NLO and NNLO corrections both vanish. For example we study tri-bimaximal (TB) mixing 
via constrained sequential dominance (CSD) which involves only a NNLO correction and 
tri-bimax;imal-reactor (TBR) mixing via partially constrained sequential dominance (PCSD) 
which involves a NLO correction suppressed by the small reactor angle and show that the 
analytic results have good agreement with the numerical results for these cases. 



* E-mail: 
^E-mail: 

* E-mail: 



antusch@mppmu.mpg.de 

sb26@soton.ac.uk 

sfk@hep.phys.soton.ac.uk 



1 Introduction 



The flavour problem of the Standard Model (SM), i.e. the origin of the observed fermion masses 
and mixings, is one of the greatest unsolved puzzles in particle physics. The discovery of neutrino 
masses and large lepton mixing angles has added another piece to this puzzle. A successful theory 
of flavour has to explain why neutrino masses are small and why the mixing in the lepton sector 
is large whereas the mixing in the quark sector is small. 

A promising approach to both issues is provided by the type I see-saw mechanism with 
right-handed neutrino dominance [2] . In particular in this paper we shall focus on the sequential 
dominance (SD) of three right-handed neutrinos [Sj SJ El El [TJ [8] . While the seesaw mechanism is 
generically realised in left-right symmetric Grand Unified Theories (GUTs) and can explain the 
smallness of neutrino masses, SD, i.e. the concept that the right-handed neutrinos contribute to 
the neutrino mass matrix via the see-saw mechanism with sequential strength, provides insight 
into how large neutrino mixing can arise in a natural way for a hierarchical neutrino mass 
spectrum {m^l > \m2\ S> \mi\ where experimentally [msl ~ 5|m2| ~ 0.05 eV 0. According 
to SD, the leading order (LO) contribution to the neutrino mass matrix comes from one single 
right- handed neutrino resulting in a single neutrino mass, namely the heaviest neutrino mass 
m3, and the "atmospheric" mixing angle 623 f2]. The second largest next-to-leading order (NLO) 
contribution to the neutrino mass matrix in SD, arising from a second right-handed neutrino, 
induces the second heaviest neutrino mass m2 as well as the "solar" and "reactor" mixing 
angles 6*12 and ^13, respectively [3]. The third largest next-to- next-to-leading order (NNLO) 
contribution to the neutrino mass matrix in SD, arising from a third right-handed neutrino is 
responsible for the lightest neutrino mass mi and this is usually considered to be sufficiently 
small so that the NNLO corrections do not perturb too much the determination of neutrino 
mixing angles in SD H El [7j . 

The analytic estimates of the mixing angles in SD have so far only been presented to LO 
O m El El C] • However the NLO corrections of order 7712/711,3 might be expected to introduce 
mass dependent corrections of order 20% to the analytic expressions for the neutrino mixing 
angles. Given the present and planned program of neutrino experiments, it is apparent that 
neutrino physics is now entering the precision era, with all mixing angles being measured in the 
forseeable future to an accuracy of a few per cent [9j. Against such experimental progress it is 
clear that theoretical errors of order 20% are no longer acceptable. This motivates a study of 
NLO corrections to SD. 

Typically unified models predict a third right-handed neutrino which contributes to the 
seesaw mechanism at NNLO and provides a mass rrii to the lightest neutrino which will also 
give some corrections to the analytic expressions for the neutrino mixing angles at NNLO. The 
importance of NNLO corrections to the neutrino mixing angles of order mi/m^ clearly depends 
on the size of mi. However, in unified models with SD, the value of mi is governed by a rather 
large third family Yukawa coupling to the third right-handed neutrino of order unity, and in such 
models the NNLO corrections could be significant. This motivates studying NNLO corrections 
to SD in addition to the NLO corrections. 

We shall refer to the NNLO corrections as "soft" since they vanish in the strongly hierarchical 
vanishing mi limit. By contrast we call the NLO corrections "hard" since it is not possible to 
decouple these corrections by taking the vanishing limit of m2. However there is a class of SD 
theories in which the NLO and NNLO corrections both vanish namely those which satisfy form 
dominance (FD) |10j (see also [H]). FD jlOj is a criterion whereby the columns of the Dirac 
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matrix are proportional to the respective columns of the neutrino mixing matrix in a basis 
where the charged lepton and heavy Majorana mass matrix are diagonal. FD implies that the 
neutrino mass matrix resulting from the type I see-saw mechanism is form diagonalizable with 
the neutrino mixing angles being independent of the neutrino masses. This observation helps to 
explain why the LO results for SD are often observed to work better than expected in particular 
cases. For example there are special cases where FD is only violated by soft corrections leading to 
NNLO corrections but not NLO corrections, namely those which satisfy constrained sequential 
dominance (CSD) [12] where a strong hierarchy is assumed and terms proportional to mi violate 



As the precision of the neutrino data has improved, it has become apparent that lepton 
mixing is consistent with the so called tri-bimaximal (TB) mixing pattern |T3] 



TB lepton mixing in particular hints at a spontaneously broken family symmetry G f which might 
underpin a flavour theory of all quarks and leptons, but which might only reveal itself in the 
neutrino sector. What is the nature of such a family symmetry? In the (diagonal) charged lepton 
mass basis, it has been shown that the neutrino mass matrix leading to TB mixing is invariant 
under the Klein group [14j. The observed neutrino flavour symmetry of the Klein group may arise 
either directly or indirectly from certain classes of discrete family symmetries [15] . Several models 
have been constructed that account for the structure of leptonic mixings, e.g. [T6 l [T7 1 ITS l [T9l |20] . 
while other models extend the underlying family symmetry to provide a description of the 
complete fermionic structure [2ll[22l[23l[a[25l[23[2Il[2Hl[Ml[30l[ai[^ 

Most of these models satisfy the conditions of either FD [10] or CSD |12j . In particular those in 
[Ml Eg EZlllHl [291 Sni EI] satisfy FD while those in[25l[32ll33llMll36ll371[38| satisfy CSD. We 
emphasise that FD and CSD are not only compatible with symmetry but are actually enforced 
by symmetry in these models. For example the texture zeros elements of the Dirac mass matrix 
that we shall encounter in Section 6 are all enforced by symmetries in realistic models. 

On the other hand there are many other models in the literature too numerous to mention 
that do not yield very precise TB mixing as a result of a symmetry for example, see |39| for 
review papers with more extensive references. Many of these other models will satisfy the looser 
conditions of SD [21 [3l [5] [U [7] but not those of CSD or FD. Such general SD models would become 
necessary if significant deviations from TB mixing are observed. For example, while TB mixing 
implies zero mixing angle ^13, the most recent experimental results hint at a non-zero (and in 
fact comparatively large) reactor mixing in the one sigma range s^g = sin^ ^13 w 0.02 it 0.01 
|40j corresponding to a reactor angle of ^13 ~ 8° it 2°. Assuming TB mixing of the neutrino 
mass matrix via CSD, various effects can generate deviations from ^13 = 0, for example charged 
lepton corrections [41j , RG running effects |421I43| H3] or corrections from canonical normalisation 
[35l H6j . However, by appealing to these effects, it is difficult to understand why the reactor 
angle should be much larger than the deviation of the solar angle from its tri-bimaximal value. 
The simplest possibility for allowing a sizeable non-zero reactor angle (e.g. ^13 ~ 8°) while 
maintaining accurately the TB predictions for the solar and atmospheric angles, 9i2 ~ 35° and 
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^23 ~ 45°, is called tri-bimaximal-reactor (TBR) mixing |47j corresponding to the mixing matrix, 
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(1.2) 



where we have introduced the reactor parameter r defined by S13 = |l8] where 3^3 ~ 0.02 
corresponds to r w 0.2. TBR mixing can arise from type I see-saw mechanism via a very simple 
modification to CSD called partially constrained sequential dominance (PCSD) ^47j. Estimates 
suggest that PCSD is capable of accommodating a sizeable reactor angle while the atmospheric 
and solar angles are predicted to remain close to their TB values [U] . However, as in the case of 
CSD, the analytic results for the mixing angles have only been presented to LO, and once again 
the NLO corrections of order m^jraj, might be expected to introduce mass dependent corrections 
of order 20%. However we shall show that for the case of PCSD such NLO corrections become 
suppressed by r leading to small corrections, for example 4% correction for r ~ 0.2. We shall 
drop such corrections in our analytic results, though they are apparent in the numerical results. 

In this paper, then, within the general framework of the type I see-saw mechanism with SD, 
we derive approximate analytic formulae of the neutrino mixing angles including both the hard 
NLO and the soft NNLO corrections. These results may be applied in a model independent 
way to a very large class of theories which satisfy SD away from the FD limit. The derivation 
of these analytic expressions builds on the results presented in |5] where the NLO and NNLO 
corrections were not considered We apply these analytic results to the cases of CSD and 
PCSD, as examples where soft and hard violations of FD are expected, respectively. However 
in PCSD the hard NLO corrections are suppressed by the small reactor angle and we drop such 
corrections in our analytic formulae. We compare the analytic results to an exact numerical 
diagonalization of the neutrino mass matrix, for the examples of two numerical GUT models 
previously studied [S] with CSD and PCSD, and obtain good agreement which support our 
conclusions based on the analytic formulae. 

The paper is organised as follows. In Section [2| we discuss lepton mixing from the special 
types of SD corresponding to CSD and PCSD as examples of soft and hard violations of FD 
leading to NNLO and NLO corrections respectively. In Section [3j we briefly review general SD 
and the LO analytic expressions for neutrino mixing angles from [5J away from the FD limit. 
Section |4] extends these results to include NLO and NNLO corrections in terms of the NLO 
parameters and the NNLO parameters In Section [5] we apply the analytic results to the 
cases of CSD and PCSD and show that the corrections correspond to soft NNLO corrections 
and hard NLO corrections suppressed by the reactor angle, respectively. In Section [6j we 
present numerical results for the mixing angles and the neutrino masses using the Mathematica 
package MPT/REAP [iS] for some examples of CSD and PCSD which support our conclusions 
based on the analytic formulae. Section [7] concludes the paper. The detailed procedure for the 
diagonalisation of the left-handed neutrino mass matrix as well as details of the derivations of 
the analytic results are given in the Appendices. 

Given the large number of acronyms introduced, it may be convenient to summarise them: 

^Although the NLO corrections were calculated for the atmospheric angle they were not considered for the 
other angles, and NNLO corrections were completely neglected [5]- 

■^Mixing Parameter Tools (MPT) is a package provided with REAP and it is mainly used to extract neutrino 
mixing parameters. 
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SM = Standard Model 

GUT = Grand Unified Theory 

LO = Leading Order 

NLO = Next-to Leading Order 

NNLO = Next-to-Next-to Leading Order 

TB = Tri-Bimaximal 

TBR = Tri-Bimaximal-Reactor 

SD = Sequential Dominance 

LSD = Light Sequential Dominance 

HSD = Heavy Sequential Dominance 

FD = Form Dominance 

CSD = Constrained Sequential Dominance 

PCSD = Partially Constrained Sequential Dominance 

REAP = Renormalisation group Evolution of Angles and Phases 

MPT = Mixing Parameter Tools 



2 Lepton mixing in special cases of sequential dominance 

The mixing matrix in the lepton sector, the MNS matrix J/mns? is defined as the matrix which 
appears in the electroweak coupling to the W bosons expressed in terms of lepton mass eigen- 
states. With the masses of charged leptons Me and neutrinos rrii, written as 

jC = -eiM^eR - ^ulitilli^I + H.c. , (2.1) 

and performing the transformation from flavour to mass basis by 

Me V^^ = diag(me, m^, m^), Kl rriLL Vj^ = diag(mi, m2, ma), (2.2) 

the MNS matrix is given by 

Umns = Kl^Jl • (2-3) 

In this paper wc shall choose a flavour basis in which the charged lepton mass matrix is 
diagonal. In this basis the MNS matrix arises from the neutrino sector, and the effective neutrino 
mass matrix rriLr, is given in terms of the neutrino masses mi, m2, ms by, 

rriLL = ^7MNsdiag(mi, m2, m3)U^^^s 

= mi$i$f + m2$2^2 + "^3^3^3 , (2.4) 

where we have written the mixing matrix in terms of three column vectors 

f/MNS = (^l,$2,^'3)- (2.5) 

Turning to the type I see-saw mechanism, the starting point is a heavy right-handed Majo- 
rana neutrino mass matrix Mur and a Dirac neutrino mass matrix (in the left-right convention) 
Mlr, with the light effective left-handed Majorana neutrino mass matrix rriLL given by the type 
I see-saw formula (up to an irrelevant minus sign which henceforth is dropped), 

ruLL = MlrM-]^MIj^. (2.6) 
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In a basis in which M^/j is diagonal, we may write, 



MRR = dmg{Y,X,X') (2.7) 
and Mm may be written in terms of three general column vectors A,B,C, 

Mlr = {A,B,C). (2.8) 

The see-saw formula then gives, 

AA^ bb^ cc^ , , 

mLL = ^ + ^ + ^. (2.9) 



By comparing Eq. (2.9) to (2.4) it is clear that the neutrino mass matrix is form diagonal- 
izable if we assume that the columns of the mixing matrix <I>i,$2,^3 are proportional to the 
columns of the Dirac mass matrix A = e<l>3, B = 6$2) C = c'^i, where e, b, d are real parame- 
ters, an assumption known as form dominance (FD) |10j (see also [H]). In this case the physical 
neutrino masses are given by 772.3 = ^ 1^ ■, ^2 = b'^ /X, mi = c'^/X' and the neutrino mass 
matrix is diagonalized precisely by Umns due to the unitarity relation ^l^j = Sij. There are no 
NLO or NNLO corrections in FD since the mixing matrix is determined by the column vectors 
$j which are independent of the parameters which determine the neutrino masses m^. This 
conclusion is independent of the choice of mixing matrix UmnSj however the usual assumption 
is that it is of the TB form. 



For example, TB mixing in Eq. (1.1) results from 





^ = 4=1, 5 = 4= i lc7 = 4= -l , (2.10) 




where we have written A = e$3^, B = b^^^ , C = c'^J^, with Utb = {^i^ ,^2^ ^^3^) 
we identify the physical neutrino mass eigenvalues as ma = e'^/Y, m2 = u /X, mi = c'^/X'. 

It is interesting to compare FD to Constrained Sequential Dominance (CSD) defined in jl2j . 
In CSD a strong hierarchy \mi\ <C \m2\ < {m^l is assumed which enables mi to be effectively 
ignored (typically this is achieved by taking the third right-handed neutrino mass X' to be very 
heavy leading to a very light mi). Thus CSD is seen to be just a special case of FD corresponding 
to a strong neutrino mass hierarchy. FD on the other hand is more general and allows any choice 
of neutrino masses including a mild hierarchy, an inverted hierarchy or a quasi-degenerate mass 
pattern. 



In practice, CSD is defined by only assuming the first and second conditions in Eq. (2.10) 
12] , with the third condition approximated as follows, 



A = 4= 1 , 5 = 4= 1 1 c«4= , (2.11 






since the Dirac neutrino mass matrix is given by M^r = {A,B,C) and in hierarchical unified 
models the third column is dominated by the large 3-3 Yukawa coupling, compared to which the 
other elements in the third column are approximately negligible. In the limit mi — )• the entire 
third column C plays no role and may be neglected, so the fact that C does not satisfy the FD 
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conditions is irrelevant, and in this limit the CSD conditions for A, B above lead precisely to the 
TB mixing angles. In particular we note that in this limit of CSD there are no NLO corrections 
to the TB neutrino mixing angles. However, in practice, the large 3-3 Yukawa coupling may be 
expected to lead to a non-zero mi, and in this case the TB mixing angles would be expected 
to be subject to NNLO corrections. CSD is therefore an example where FD is only violated by 
soft NNLO corrections. The neutrino mass mi can be written approximately at CSD, in terms 
of the 3-3 Yukawa coupling, as |^ 

™. - (2.12) 



Assuming a strong neutrino mass hierarchy |mi| <C |m2| < |m3|, TBR mixing in Eq. (1.2) 
results from a simple modification to CSD, corresponding to allowing a non-zero 1-1 element of 
the Dirac neutrino mass matrix. 





A = —\ 1 \, B = —\ 1,C«— , (2.13) 



where we have written A = e^s, where <^3 is the third column of the TBR matrix in Eq. (1.2). 
This is referred to as Partially Constrained Sequential Dominance (PCSD) |37|, since one of 
the conditions of CSD is maintained, while another one is violated. We emphasise that B is 
unchanged from the case of CSD, and in particular B is not proportional to the second column 



of the TBR matrix in Eq. (1.2), and thus, as noted in [37] . PCSD violates FD since A and B are 



not orthogonal if ^13 is non-zero. Thus, one might expect that the TBR form of mixing matrix 



in Eq. (1.2) will result only to LO with NLO corrections of order |m2|/|m3| generally expected. 
However the unitarity of the columns A, B which, if exact, would imply no NLO corrections, is 
only spoiled by the first element of A proportional to the reactor angle. Therefore we expect 
the NLO corrections to PCSD to be suppressed by ^13, and our numerical results confirm this. 

We note that the special cases of SD considered here, namely CSD and PCSD, have a 
certain theoretical elegance since the columns of the Dirac neutrino mass matrix A, B take 
very simple forms in these cases. As discussed in the Introduction these simple forms may be 
achieved in flavour models based on non-Abelian discrete family symmetry (containing triplet 
representations) and vacuum alignment of the flavons which break the family symmetry. In 
particular the CSD and PCSD approaches considered here correspond to flavour models in 
which the neutrino flavour symmetry inherent in TB and TBR mixing is achieved in an indirect 
way [15] from the family symmetry as an accidental symmetry. In such indirect models [25\ [32l 
\33\ [Ml l36l [371 138] the flavon vacuum alignment responsible for the simple forms of the columns 
A,B of the Dirac neutrino mass matrix arises in supersymmetric models from the so called 
D-term vacuum alignment mechanism |50j . 

In general, the MNS matrix may not be as simple as either the TB or the TBR patterns 
suggest and thus the special cases of SD such as CSD or PCSD may be regarded as being too 
restrictive. Indeed it is quite possible that the actual lepton mixing angles deviate significantly 
from those predicted by these special cases, and thus the hierarchical neutrinos satisfy more 
general SD conditions (specified more precisely in the next section) than the CSD or PCSD 
ones. In such cases it is useful to have analytic expressions for the neutrino mixing angles for 
the general case of SD, including the NLO and NNLO corrections, due to arbitrary violations of 



^Another formula for the mass mi is derived in Appendix [b1 
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FD in which the columns A,B,C do not satisfy complex orthogonality, not even approximately, 
and these will be given in subsequent sections. 



3 Neutrino mixing angles in general SD to LO 

We consider the case of see-saw mechanism with general SD, involving a right-handed neutrino 
mass matrix M/j/j and a Dirac neutrino matrix M^r. We take Mrr to have a diagonal form 
with real eigenvalues as follows. 



Mrr =0X0. (3.1) 




We also write the complex Dirac neutrino mass matrix M^r in terms of the general (uncon- 
strained) mass matrix elements a, b, c, d, e, /, a', b' , c' as, 

(d a a' \ 
e b b' \. (3.2) 
f c c' J 

In this paper, motivated by unified models where the 3-3 element of the hierarchical mass matrix 
is very large, we take the two complex elements a', b' to be negligible in comparison to c'. The 
effective neutrino mass matrix can be derived using the see-saw formula [1], 

AA^ BB'^ CC'^ 
rriLL = MlbMrrMIr = ^ + + ^ ' (3-3) 

which gives the following complex symmetric matrix, 



■mLL 



/ a?_ j_ iP_ ab , de ae , df 

I X ^ Y X ^ Y x"|"y 

ab I de fe^ I be i ef I ^o A\ 

X + Y X+Y X+Y \ (3.4) 

ae , df be , ef c'2 , , 



X ^ Y X' ^ X ^ Y 

This mass matrix can be diagonalised, for the case of hierarchical neutrino mass matrix, as 
presented in Appendix [A} 

Sequential dominance corresponds to the following condition, ignoring d to start with, 



\ms\ = O ^' J » |m2| = O ^-^j » |mi| = O (^^J (3.5) 

where x,y £ a,b, c and all Dirac mass matrix elements are assumed to be complex. Diagonalising 

the neutrino mass matrix, the mixing angles can be derived to leading order (LO), in the 
framework of the seesaw mechanism, as presented in [5l |6] 

t23 ~ ^ = ^23, (3.6) 

~ — TTi /,/J ^ n 7jJ\ = ^12' (3-7) 

C23|0| COS(0'J - S23|c| COs((/)'c) 



^ \a\ie*b + rc) y _ 



(3.8) 



Throughout this paper, for brevity we often write Cij = cos(^jj), Sij = siii(^ij), tij = taii(^ij) 
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where some of the Dirac masses were written 

|g«<Pa:_ 'pj-^g phases x ^-nd (j)2 are fixed, to 

give real 612 and ^13 angles, by: 



where 



C23\b\sm{(j)'f,) ^ S23|c| sin((/)'c), 

(j)2 - (t)a - 4>* ■, 

(t>b = <t>h - (t>a - (t>2 - X, 

(p* = arg{e*b + f*c). 
In the large d limit, the angle ^13 can be expressed as follows [5]: 

\d\ 



713 



V\e\' + \f\' 



90 

^13- 



(3.9) 
(3.10) 



(3.11) 
(3.12) 
(3.13) 



(3.14) 



Note that ^13 and ^5*3 are given differently in the small d and large d cases so we must be careful 



to distinguish the two limiting cases. The phases 02 and (^3 appearing in Eq. (A. 9) are fixed by: 



4>2 = 4>e- (Pd 

(p3 = (Pf - 4>d- 



(3.15) 

(3.16) 



As a simple application of these results, it is clear that TB neutrino mixing (tan023 = 1; 
tan 012 = l/\/2 and ^13 = 0) can be achieved by considering CSD which corresponds to the 
following set of conditions for the Yukawa couplings |12j . 



a 


= \b\ = \cl 


(3.17) 


\d\ 


= 0, 


(3.18) 


e 


= I/I, 


(3.19) 


e*b + f*c 


= 0. 


(3.20) 



These conditions are also in accordance with Eq. (2.11). 

The results in this section were derived in [5j £j at LO, neglecting the NLO and NNLO 
corrections. For the remainder of the paper, we will refer to the LO mixing angles as ^23)^12 ^^'^ 

respectively, as indicated in Eqs. (3.6), (3.7), (3.8) or (3.14). Note that the LO atmospheric 
and solar neutrino mixing angles do not depend on the neutrino masses, while the reactor angle 
(in the small d case) does. However there are expected to be mass dependent corrections to all 
angles at order NLO and NNLO. As mentioned, near to the CSD limit, the NLO corrections 
cancel, leaving only the NNLO corrections arising from the 3-3 Dirac neutrino mass matrix 
element c'. This implies that, providing that the NNLO corrections are not too large (i.e. mi 
is not too large) the simple LO approximations to the neutrino mixing angles presented in 
this section are more accurate than might be expected in the case of CSD. More generally the 
LO results may be sufficient for examples of SD close to the FD limit where NLO and NNLO 
corrections vanish. However, for many other cases of SD away from the FD limit, the LO results 
are not sufficient and the NLO and NNLO corrections need to be considered. This is done in 
the next section. 
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4 Neutrino mixing angles in general SD to NLO and NNLO 



In this section, we derive approximate analytic expressions for neutrino mixing angles in the case 
of neutrino mass hierarchy in general SD including NLO and NNLO corections. The derivations 
make use of the diagonalisation procedure outlined in Appendix Note that this procedure 
relies on the fact that the 1-3 reactor angle is small which enables the diagonalization to be 
performed in three stages: first diagonalize the 2-3 block, then the 1-3 outer block, and finally 
the 1-2 block. Such a method cannot take into account subleading corrections resulting from 
changes to the 2-3 mixing angle as a result of 1-3 rotations and therefore NLO corrections to 
the 2-3 angle suppressed by the 1-3 angle are not present in the analytic results. 



4.1 The atmospheric angle 

As discussed in Appendix |Aj the diagonalisation of the mass matrix involves applying the real 
rotation R23 after re-phasing the matrix. This rotation gives rise to two new mass terms 77122 



m'^ given by Eqs. (B.5),(B.3) respectively. We consider the 23 block as presented in Eq. (A. 10) 
which gives rise to an expression for tan(2023) in terms of the lower block masses 771-23, ^-22 and 



77733. Using Eq. (3.4), we can substitute for these masses in terms of the Yukawa couplings to 



get the following expression, 
tan(2023 



2^(l + ei)e^(- 



■<t>2-4>3) 



f (1 + 62 + r/i)e^(-203) -^{1 + e3)e^(-202) ' 
where we have introduced new parameters ei, £2, £3 and rji, which are given as follows, 



(4.1) 



be 
X_ 

£/ 
Y 



£2 



c 

X 

/!' 

Y 



€3 



Y 



c 

X' 

11 
Y 



(4.2) 



Note that ei,rji are of order 7^2/7773,7771/7773 respectively, so that parametrize the NLO cor- 
rections while rji parametrize the NNLO corrections. 

Introducing the small parameter 5 such that |/| = |e|(l — 5), we get our final result for the 
atmospheric angle in SD: 

i23«4(l+ Re(7)), (4.3) 
where the complex couplings e, / are written in terms of their absolute values and phases as 

complex parameter 7 is written as: 



e = |e|e*'^% / = |/|e*'^-f, respectively, t23 = tan(023)|e,=o,7yi=o is given by Eq. (|3.6|) and the 



1 , 5 . , 

7 ~ 2*-^^ - €2 - Vl) + 2(^3 + £2 - 2ei + rji). 



(4.4) 



4.2 The reactor angle 

We apply the i?i3 rotation, as outlined in the Appendix, in order to modify the outer block of 
the mass matrix. We consider the reduced matrix that only involves the 13 elements and this 



gives rise to two zeros in the 13, 31 positions as presented in Eq. ( A.ll ). The neutrino angle ^13 
can then be written as, 



1 



^13 



771 



/O 



K3(l-7(s' 



^3)^) 



+ e" 



-l<p2 



ah de..,^ ^. 



(4.5) 
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where the masses mf,mi^, are given by Eqs. (B.4), (B.IO) respectively. The complex parameter 
P is given by: 

/3 « i43?e3 + {4s)\e2 + vi) - e,e-^'^% (4.6) 
where the NLO correction parameter £4 is defined as, 



64 



X 



Y 



We can simplify Eq. (4.5) further, after expressing the masses in terms of the complex couplings, 



by considering two different limits, namelyt the large d limit and the small d limit, as follows: 
• In the large d limit, ^ >> the angle ^13 can be expressed as, 

ei3«0?3(l- Re (/?)), (4.7) 



where the angle 05*3 = ^i3|r?j=o,ei=o is given by Eq. (3.14) 



In the small d limit, i^?^, ^ << which is usually the case in CSD, ^13 can be 



expressed as. 



?13 « 0?3(1- Re(7)(s 



^2 
23 J 



Re(/3)) 



+ ■S23k5l ((Re7cos((/>') - Im7sin(0'))^ + (Re7sin((/)') + Im7cos(</>'))^) ^ (4. 



where 4>' = (j)2 — "^(pe and 9^^^ in this limit, is derived in ^ and given by Eq. (3.8) and the 
NLO correction parameter £5 is defined as. 



£5 



ah (\e\^ + \fr-' 



X 



Y 



(4.9) 



In the PCSD case with non-zero d, we can write the leading result for ^13 as, 



713 



??3 + 



\d\ 



[\ - Re (7)(4)' - Re (/3)) + 4 Re (7). 



|d||e| 



+ ■523|e5| ((Re7cos((7:)') - Im7 sin((/.'))^ + (Re7sin(</.') + Im7cos((/)'))^) ^ ,(4.10) 



where and the parameter €5 are given by Eqs. (3.8), (4.9) respectively. 
4.3 The solar angle 



As shown in Eq. (A.6), applying the phase matrix Pi introduces a new phase x to the mass 



matrix. We can then apply the real rotation i?i2, as presented in Eq. (A. 12), which modifies the 



matrix by putting zeros in the 12,21 positions. Using Eqs. (B.5), (B.7) and (B.ll), we get the 
following expression for tan(26'i2), 



tan(26' 



12J 



2AB 



(l-7(sy'-Ci-C2(l-Ci-7(s^3)')), 



(4.11) 



11 



where, similarly to [5], A, B are expressed in terms of the complex Yukawa couplings as, 

A = ^, (4.12) 



a 
X 



, h — rp'^i<t>e-<Pf) 



X 



(4.13) 



and the new parameters C,i and C,2 are given, in the small d limit, to first order in 7 and /3 as, 



Ci 

C2 



p-i{<t>i+x) ( "^'^23 A 

\ABX)^ 

U2 „2 



1 



52 -yl2 



??2 



where 772 is given by. 



.J2 /|,|2 + |^|2x-l 



Using Eq. (4.11) we obtain our final result for the solar mixing angle in SD: 
where = tan(0i2)|,;j=o,ei=o is given by Eq. (3.7). The new parameter C,' is given as, 



(7(4)' + Ci + C2(i-Ci-7(4)'))- 



(4.14) 
(4.15) 

(4.16) 

(4.17) 

(4.18) 



52 + ^2 

5 Neutrino mixing angles in special cases of SD 

In this section, we will look at how the rather complicated analytic results for the mixing angles 
in SD, derived in the previous section to NLO and NNLO, can be simplified in the special cases of 
SD corresponding to CSD and PCSD. For simplicity, we shall take (/>e = tt and all the remaining 
Yukawa phases are taken to be zero except 0^/ which is left general. 

5.1 CSD 



CSD corresponds to SD with the constraints defined in Eq. (3.17). A particular example of 



CSD was discussed in Eq. (2.11) below which we presented an argument for why we expect 



the NLO corrections to vanish for CSD leaving only the NNLO corrections dependent on mi 
whose magnitude is governed by the 3-3 element of the Dirac neutrino mass matrix c' . Using the 
analytic results for neutrino mixing angles in SD derived to NLO and NNLO we can explicitly 
verify that the NLO corrections vanish for CSD. 

5.1.1 The atmospheric angle 



The atmospheric angle in Eq. (4.3) becomes, in CSD. 



23 



1 + Re (7 



CSD\ 



(5.1) 
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which involves a correction 7 in Eq. (4.4) which depends on the NLO parameters and the 



NNLO parameters r]i in Eq. (4.2). In CSD the conditions in Eq. (3.17) imply that the are all 
equal and 6 = 0. From Eq. (4.4) it is clear that the NLO contributions to 7 described by the 
€i cancel. This implies that the atmospheric angle is corrected by 7 which only involves NNLO 
corrections given by, 



CSD 



m 

' 2 



1 |cf F 
'2 \e?X' 



(5.2) 



7 

5.1.2 The reactor angle 

Turning to the reactor angle ^13, we only need to consider the expression valid in the small d 
limit given by Eq. (4.8). Imposing the CSD conditions in Eq. (3.17), the angle ^5*3 becomes 



exactly zero as can seen from Eq. (3.8). As a result, the first term of Eq. (4.8) vanishes. Also 



the third term vanishes for CSD. We are only left with the second term of order r]e, 



qCSD 
^13 



cos(2(/)(. 



(5.3) 



' 4^/2 WXX' 

This implies that the reactor angle is given by a term proportional to NNLO • NLO corrections. 



5.1.3 The solar angle 



The solar angle in Eq. (4.17) becomes, in CSD 

1 



i-CSD 
''12 



^/2 



(l-Re(C'^^^)), 



which involves a correction ^' in Eq. ( 4.18| ) which we approximate here to, 

iCSD !i ( J |_ ^CSD j_ /-CSD \ 



: ^CSD . /-CSD 

^ + Ci +C2 J- 



(5.4) 



(5.5) 



which depends on ^^^^ as well as the parameters Ci, C2 in Eqs. (4.14), (4.15), which also depend 



on /3 in Eq. (4.6). The parameter ^^^^ takes the simplified form above in Eq. (5.2). The 



parameters /3,Ci,C2i given by Eqs. (4.6), (4.14), (4.15), can also be simplified in CSD as 



I^CSD ^ 




/-CSD ^ 


^CSD 

2 


/■CSD ^ 
<,2 


l|cfX 
2\a\^X' 



(5.6) 
(5.7) 

(5.8) 



The important point to note is that the solar angle is corrected only by NNLO corrections. 
Indeed all corrections to neutrino mixing angles in CSD vanish at NLO, and first occur only at 
NNLO, as anticipated. 

We have compared the analytical results in CSD with the numerical results calculated using 
the MPT/REAP software package [IH] and found good agreement. The numerical results are 
shown in Fig. [l] (left panel). While ^13 is rather stable and remains close to zero, the plot 
in Fig. [T] illustrates that the deviation of 9i2 from the tri-bimaximal value is about twice the 
deviation of ^23 from 7r/4. It also illustrates the above-derived dependence of the results on the 
phase (pel- 
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5.2 PCSD 



PCSD is similar to CSD defined in Eq. 
example of PCSD was discussed in Eq. 



3.17), but with a non-zero value of d. A particular 
. In this example the NLO corrections do not 
e. However, in deriving the analytic results for 



2.13 



vanish but are suppressed by the reactor angL 
neutrino mixing angles in SD to NLO and NNLO, such corrections are beyond the scope of 
our approximations. However the results below have the simple property that, as in the case 
of CSD, they only depend on NNLO corrections. This demonstrates that the NLO corrections 
vanish for PCSD, neglecting NLO terms suppressed by the reactor angle. 

5.2.1 The atmospheric angle 



The atmospheric angle in Eq. (4.3) becomes, in PCSD, 

t^3'^5^;^l+Re(7^^^^), (5.9) 

where it is easy to see that the result is identical to the case of CSD, 

^PCSD ^ ^CSD_ (5_10) 

This implies that the atmospheric angle correction only involves NNLO corrections, as in the 
case of CSD. 

5.2.2 The reactor angle 

Turning to the reactor angle ^13, we find for PCSD, using Eq. (|4.5l). 



^fs^^^ - ^?3(1 + Re (7^^^^)) - ''"K''''^ ^Zx ' ^'-''^ 



where the LO expression is given by the large d approximation in Eq. (3.14) 



no ^ \d\ 



(5.12) 



The reactor angle therefore only receives NNLO corrections proportional to ^^'^^^ = r^csD ^ 
5.2.3 The solar angle 



The solar angle in Eq. (4.17) becomes, in PCSD, 



tf2^"^-^(l-Re(r^^^)), (5.13) 



which involves a correction C,' in Eq. ( 4.18[ ) which we approximate here to 



PCSD 

(5.14) 



^iPCSD ^ ^ J_ L /-PCSD . /-PCSD \ 



14 



50 


023 


1 




50 


623 ; 


40 


■ 0,2 




<t>d=0 ~ : 


40 


■ 0,2 - 


30 


- 






30 




20 








20 




10 








10 







0-\3 




















-4.6 -4.4 -4.2 -4.0 -3.8 -3.6 -3.4 o.OO 0.05 0.10 0.15 0.20 0.25 0.30 

log,„(m,,,/eV) |d|/|e| 



Figure 1: Numerical results using the MPT/REAP package |49] for the corrections to the lepton mixing angles in 
CSD (left panel) and PCSD (right panel), where in the latter case c' has been set to zero. In CSD (left panel), the 
NNLO correction depends on the mass of the lightest neutrino, m^j , and the light-grey vertical band shows the 
typical expectation for m^^ in some classes of GUT models. While ^13 is rather stable and remains close to zero, 
the deviation of Q\2 from the tri-bimaximal value is about twice the deviation of ^23 from 7r/4. The two extreme 
cases with = 0, 7r/2 (note that the change of both, Q\2 and S23 is approximately proportional to cos(2<j!>c/)) are 
displayed and the grey area between the curves corresponds to other values of <j!>c'. In PCSD (right panel) with 
zero c', Q\2 and 6123 are nearly unchanged and only &13 increases with increasing |rf|/|e| = |d|/|/|. The results 
agree well with the analytic results, and confirm that in the case of PCSD there there are small NLO corrections 
to 6^23 and 6112 proportional to ^13 which are not seen in the analytic results. 



where 



^PCSB ^ ^^^^^(1 + ^03^^!^), (5.15) 
Cr^^ « C2^^^ + ^/2e?37^^^^ + (e?3)^ (5.16) 



and where 05*3 is given in Eq. (5.12, The important point to note is that the solar angle has no 
NLO corrections. However there is a correction of order (^13)^. This can be understood if one 
thinks of the diagonalization procedure since the 2-3 diagonalization in PCSD is the same as in 
CSD (since the 2-3 blocks are identical), while the subsequent 1-3 diagonalizaton will correct 
the 1-1 element at order (^13)^, and this will lead to the above correction to the solar angle at 
order (^13)^- However the results show that the corrections to neutrino mixing angles in PCSD 
vanish at NLO, neglecting NLO terms suppressed by the reactor angle, with only the NNLO 
corrections remaining to this approximation. 

The analytical results in PCSD agree well with the numerical results calculated using the 
MPT/REAP software package |l9]. For the special case of zero d the numerical results are 
displayed in Fig. [T] (right panel). In accordance with the analytical results, Q\2 and ^23 are 
nearly unchanged and only ^13 increases with increasing |(i|/|e| = 
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6 Analytical vs. numerical results in two example models 



In the previous section, analytic expressions of the mixing angles, involving NLO and NNLO 
corrections, were derived. Approximate results in the case of CSD and PCSD were also presented 
and the NLO corrections vanished in both cases. In this section, we evaluate the analytic results 
for two different numerical GUT inspired models. These are models of light sequential dominance 
(LSD) and heavy sequential dominance (HSD) [21 17] previously studied in [44J. We also present 
numerical results for the neutrino masses, presented in Appendix [B} as well as the difference in 
neutrino masses squared for both models. 



6.1 Results for the LSD model 
6.1.1 The CSD case 

We consider a simple numerical example with a diagonal right-handed neutrino Majorana mass 
matrix M^j/j given by, 

/ 5.1 X 10-^ \ 

Mrr = 7.05 X 10-9 M3, 

V 01/ 

where M3 = lO^^Ge^. We also consider a neutrino Yukawa matrix of the form, 



(6.1) 




(6.2) 



where we have taken the complex coupling d to be zero. As required by CSD, we take \a\ = 
\b\ = \c\ = 8.125 X 10-5 and |e| = |/| = 2.125 x 10"^. The value of the third family coupling c' 
is taken to be |c'| = 0.5809. We choose all the phases of the Yukawa couplings to be zero except 

(j)e {(j)e = Tt). 



Parameter 


\d\ 


|c'| 


^23 (°) 


^13(°) 


0i2n 


nil (eV) 


7712 (eV) 


7713 (eV) 


Analytic 





0.5809 


44.44 


0.04 


33.75 


1.52 X 10"^ 


8.8 X 10-3 


0.055 


MPT/REAP 





0.5809 


44.38 


0.05 


33.69 


1.61 X 10"* 


8.8 X 10-3 


0.054 



Table 1: Analytical vs. numerical results for the mixing angles and masses, evaluated in the CSD case with c' 7^ 0, 
for a model with light sequential dominance. 

Numerical results for the mixing angles, evaluated in CSD using the analytic formulae, are 
presented in Tab.[T] This table also shows numerical results obtained using MPT/REAP package 
|39] . which appear to be very close to the ones obtained through the analytic approach. We note 
that here and in the remainder of the paper, the MPT/REAP results were evaluated using the 
MPT package without considering RG running. As can be seen from Tab.[l} all the values of the 
mixing angles are slightly deviated from their TB values and this is mainly due to the presence of 
the non-zero 3-3 Yukawa coupling c'j^ In addition, we present numerical results for the neutrino 
masses mi,m2 and given by Eqs. ( |B.13 ), (B.14), (B.15), using both MPT/REAP and the 
analytic formulae. As presented in Tab. [T| we can see that the MPT/REAP and the analytic 
results are very close particularly in the case of m2- 



'In the limit c' = 0, the analytic results give exact TB values (^23 = 45.00°, 612 = 35.26°, 61 



0.00° 
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6.1.2 The PCSD case 



We consider the previous LSD model in the case of PCSD with non-zero Yukawa coupling 
d = 0.2|e|, |e| = 2.125 x 10~^ and \c'\ = 0. Keeping all the other conditions of CSD satisfied 



as outlined in Section 6.1.1, we found that the numerical values of all the mixing angles are 
deviated from their TB values particularly the reactor angle ^13 which becomes larger than zero 
and takes a value of 8.22° as shown in Tab. [2j This large value satisfies the predictions of TBR 
mixing and it is in good agreement with the most recent experimental results [;40j. 

MPT/REAP results for the neutrino mixing angles in this case are slightly different than 
the analytic results as presented in Tab.[2j This is mainly due to the approximate nature of the 
diagonalisation procedure that we followed in this paper. Tab. [2] also shows numerical results 
for the neutrino masses mi, 7712 and evaluated in PCSD using both MPT/REAP and the 
analytic expressions. As expected, the neutrino mass mi is exactly zero in this case due to the 
vanishing NNLO corrections. We note that the results for the masses m2,m^ in the analytic 
case are slightly different than the MPT/REAP case. This is also due to the difference in the 
diagonalisation approaches considered in the two cases. 



Parameter 


\d\ 


|c'| 


^23(°) 


%(°) 


^12(°) 


mi (eV) 


m2 (eV) 


ms (eV) 


Analytic 


0.2|e| 





45.00 


8.10 


35.08 





8.5 X 10"^ 


5.38 X 10"^ 


MPT/REAP 


0.2|e| 





44.29 


8.53 


34.89 





8.4 X 10"^ 


5.4 X 10"^ 



Table 2: Analytical vs. numerical results for the neutrino mixing angles and masses, evaluated in the PCSD case 
for a model with light sequential dominance, with c' = Q and \d\ — 0.2|e|. 



In order to compare our numerical values to experimental data, we present numerical results 
for the difference in the squares of neutrino masses Am^^^ and Am^^^, evaluated for the LSD 
model in CSD and PCSD, as shown in Tab.[3j These results are evaluated at the SD cases using 
both the analytic results as well as MPT/REAP. The numerical results, as shown in Tab. [3| are 
within the most recent experimental ranges presented in [30] particularly for the value of Am^^^ 
in CSD which is very close to the best fit value of 7.6 x lO~^eV^. 



Parameter 


Analytic 


MPT/REAP 


Analytic 


MPT/REAP 


SD hmit 


CSD 


CSD 


PCSD 


PCSD 




7.5 X 10"^ 


7.5 X 10"^ 


7.3 X 10"^ 


7.1 X 10-^ 




2.11 X 10-3 


2.04 X 10-3 


2.05 X 10-3 


2.1 X 10-3 



Table 3: Analytical vs. numerical results of the difference in the squares of neutrino masses (Am^oj and Am^^^) 
evaluated for the LSD model. The results are presented at CSD with non-zero c' as well as the PCSD case with 
zero c' and non- zero coupling \d\ — 0.2|e|. 



6.2 Results for the HSD model 

To check the generality of our numerical results, we consider another model with heavy sequential 
dominance (HSD). The right-handed neutrino Majorana mass matrix Mrr, in this case, is given 
by, 

/ 3.991 X 10^ \ 

Mrr = 5.8 X 10^° . (6.3) 

V 5.021 X W'^ I 
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This model satisfies HSD where the dominant contribution to the neutrino mass is coming 
from the heaviest right-handed neutrino. The neutrino Yukawa matrix is of the form given in 
Eq. (6.2) with the following values of the Yukawa couplings: \a\ = \b\ = \c\ = 2.401 x 10~^ 

-5 



a\ = \b\ = \c\ 
Similarly to the LSD model, we take all the phases of 
vr). 



|e| = I/I = 0.677 and |c'| = 2.992 x 10 
the Yukawa couplings to be zero except the coupling e (0, 

Analytic and MPT/REAP results of the mixing angles and masses, in CSD, are presented 
and compared as shown in Tab. [4} We note that, for this model, the values of the mixing angles 
are much closer to their TB values compared to the LSD model, which is mainly due to the 
smallness of the 3-3 Yukawa coupling c' in this case. We also present results for the PCSD case 
with non-zero 1-1 Yukawa coupling (d), as shown in Tab. [sj Analogously to the LSD model, 
the reactor angle for the HSD model, is found to be large and within the recent experimental 
range presented in [30]. The neutrino mass mi is exactly zero at the PCSD case with c' = as 
expected. 



Parameter 


\d\ 


W\ 


^23 (°) 


^13(°) 


0i2n 


mi (eV) 


1712 (eV) 


?Ti3 (eV) 


Analytic 





2.99 X 10-5 


44.96 


0.003 


35.18 


1.01 X 10-^ 


9 X 10-3 


0.055 


MPT/REAP 





2.99 X 10-^ 


44.96 


0.003 


35.16 


1.1 X 10-^ 


9 X 10-3 


0.055 



Table 4: Analytical vs. numerical results for the neutrino mixing angles and masses, evaluated in CSD with 
c' 7^ 0, for a model with heavy sequential dominance. 



Parameter 


\d\ 


|c'| 


^23(°) 


^13(°) 


^12(°) 


mi (eV) 


m2 (eV) 


ma (eV) 


Analytic 


0.2|e| 





45.00 


8.10 


35.08 





9 X 10-3 


0.055 


MPT/REAP 


0.2|e| 





44.27 


8.55 


34.89 





8.9 X 10-3 


0.056 



Table 5: Analytical vs. numerical results for the neutrino mixing angles and masses, evaluated in the PCSD case 
for a model with heavy sequential dominance, with c' = and \d\ = 0.2|e|. 



Similarly to the LSD model, we present numerical results for the difference in the squares 
of neutrino masses Am^^^ and Am^j^, evaluated for the HSD model, as shown in Tab.|6j The 
results for this model, which are also presented at both SD cases using analytic results as well 
as MPT/REAP, are within the most recent experimental ranges presented in |40j . We note that 
the values of Am^^^, in all cases presented in Tab. [6| are closer to the upper limit of the 3a 
experimental range [40J. 



Parameter 


Analytic 


MPT/REAP 


Analytic 


MPT/REAP 


SD hmit 


CSD 


CSD 


PCSD 


PCSD 


Am^ (eV) 


8.2 X 10-5 


8.15 X 10-5 


8.2 X 10-5 


8 X 10-5 




2.16 X 10-3 


2.13 X 10-3 


2.16 X 10-3 


2.2 X 10-3 



Table 6: Analytical vs. numerical results of the difference in the squares of neutrino masses (Am^oj and Am^j„) 
evaluated for the HSD model. The results are presented at CSD with non-zero c' as well as the PCSD case with 
zero c' and non- zero coupling |d| = 0.2|ej. 

To summarize, the comparison between the analytical results and the numerical values using 
MPT/REAP has shown small differences, which are however within the expected range due to 
the approximate nature of the diagonalisation procedure followed in this paper. Explicitly, 6*12 
from the analytical results has been larger by about 0.2° than the MPT/REAP value, ^23 is 
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larger by 0.7° while the reactor angle ^13 is smaller by about 0.4° than the MPT/REAP value 
for both models. We have also presented numerical results for the neutrino masses as well as 
the difference in the squares of neutrino masses (Am^^^ and Am^^^), in the SD cases, using 
MPT/REAP and the analytic formulae. The numerical results for Am^^^ and Am^^^ are within 
the most recent experimental ranges |40j . 

7 Summary and Conclusions 

SD provides a natural way of having large mixing angles together with hierarchical neutrino 
masses from the type I see-saw mechanism. In this paper we have derived analytic expressions 
for the neutrino mixing angles including the NLO and NNLO corrections arising from the sec- 
ond lightest and lightest neutrino masses. We have pointed out that for special cases of SD 
corresponding to form dominance the NLO and NNLO corrections both vanish. This observa- 
tion helps to explain why the LO results for SD are sometimes observed to work better than 
expected in particular cases. For example, we have studied tri-bimaximal mixing via constrained 
sequential dominance which involves only a NNLO correction as an example of softly broken 
FD where the corrections vanish in the strong hierarchical limit. We have also considered tri- 
bimaximal-reactor mixing via partially constrained sequential dominance and shown that this 
involves NLO corrections to neutrino mixing angles suppressed by the small reactor angle which 
renders such corrections negligible. This supports the PCSD proposal as a robust way of obtain- 
ing a significant non-zero reactor angle with the atmospheric and solar angles remaining close 
to their tri-bimaximal values to an accuracy of a few per cent. 

We have evaluated these analytic results for two explicit GUT inspired models of so-called 
LSD type and HSD type including, in the case of CSD, the deviations from exact TB mixing 
due to the presence of a non-zero third family Yukawa coupling. For both models the analytic 
results agree well with the numerical results obtained from extracting the mixing angles using 
the MPT tool provided with MPT/REAP. In the CSD case the absence of NLO corrections is 
confirmed numerically while the soft NNLO corrections are observed. We have also seen that 
the analytic expressions in the case of PCSD with zero third family Yukawa coupling involves 
no NLO corrections in the approximation that NLO corrections multiplied by the reactor angle 
are dropped. The numerical results for PCSD confirm that the atmospheric and solar angles 
remain very close to their TB values while the reactor angle, ^13, is much larger than zero, and 
confirm that the NLO corrections are suppressed by the reactor angle. If a sizeable reactor 
angle is measured it may therefore be possible to retain the tri-bimaximal predictions for the 
atmospheric and solar angles to a few per cent accuracy using PCSD. We note, however, that 
this neglects other model dependent charged lepton, RG and canonical normalization corrections 
which have been considered elsewhere. 

To conclude, we have derived approximate analytic results which allow to extract the neutrino 
mixing angles in the framework of see-saw mechanism with sequential dominance including NLO 
and NNLO corrections from m2 and mi. We have shown that, for the special case of FD, the 
NLO and NNLO corrections vanish. CSD is an example where FD is only violated softly leading 
to NNLO corrections only, a result supported by numerical evaluation. Such NNLO corrections 
are typical in many classes of GUT models with hierarchical Yukawa matrices which lead to LSD 
with mi arising from the large third family Yukawa coupling. We emphasise that in realistic 
models, FD and CSD are enforced by family symmetry and so in testing these schemes one is 
really testing particular models based on family symmetry. It is possible for FD to yield either 
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a normal mass ordering or an inverted one, while CSD always corresponds to a normal mass 
hierarchy. We have also shown that, in the case of PCSD, the NLO corrections are suppressed 
by the reactor angle, and have again verified this result numerically. In other special cases of 
SD which are close to the FD limit, one can similarly expect the NLO and NNLO corrections 
to neutrino mixing angles to be suppressed. However, if significant deviations from TB mixing 
are observed, then the more general cases of SD away from the FD limit, as considered in this 
paper, would become relevant. The analytic results for neutrino mixing angles in SD presented 
here, including the NLO and NNLO corrections, are therefore applicable to a wider class of 
models and may be reliably used when confronting such models with high precision neutrino 
data, providing useful insights beyond a purely numerical analysis. 
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A Diagonalisation of left-handed neutrino matrix 

In this Appendix, we will briefly review the procedure of diagonalising the neutrino mass matrix 
following [5] closely. We start by writing the left-handed neutrino mass matrix as, 

77711 ini2 
^LL = I "^12 "T-22 
77713 "T-23 

In general, we diagonalise a complex, hierarchical, neutrino matrix by following a sequence of 
transformations [5], 

/ 7771 \ 

PrRi2^PrKfR2fP2'*rnlLP2'R2lRilPi'Ri'2P3' = 7n2 , (A.2) 

V 7773 / 

where the resulting matrix includes the three different neutrino masses 7n,i, 7n,2 and 7713. Rij, i,j = 
{1,2,3} are a set of real rotations, involving the Euler angles Oij, which can be written as, 

(A.3) 
(A.4) 
(A.5) 
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The matrices Pi in Eq. (A. 2) are the phase matrices, involving the phases (j)2, (ps, X Ui, which 
we write as, 



1 














1 


1 




e'^' 







e 







e 











(A.6) 

P2 = I e^*^ I (A.7) 

/ e*^i \ 
P3 = e*^2 (A.8) 
\ e*^3 / 

We briefly summarise the different steps of diagonalisation following [5]. We start by multi- 



plying the mass matrix, given by Eq. ( A.l ), by the inner complex phase matrices P2, which helps 
in performing the diagonalisation procedure. This process modifies the phases of the matrix as 
follows, 

Imiile^"^" |mi2|e*(<^i2-</'2) |mi3|e^(*i3-03) 

-<^2) |^22|e*(<^22-202) |m23|e*('^23-02-03) I (A. 9) 



friLL = \mi2\e'- 



|^^3|e»Wl3-</'3) |^23|e*(<^23-<^2-03) |j„33|ei{033-203) 



After re-phasing the matrix, we proceed by applying the real rotation i?23, defined in Eq. (A.3). 
This step modifies the lower 23 block of the mass matrix by putting zeroes in the 23, 32 elements 
of the matrix 

^22 \ _ ( |m22|e^('^22-2<^2) |^23|e*('^23-<^2-03) \ 



m', j - ^23 1^ |^^^|gife-<A2-<^3) |m33|e*('^33-203) ) ^23 l^-^Uj 

This diagonalisation not only modifies the masses m22 and 77133 but also all the other mass 
entries except mn. Substituting for these masses in terms of the Yukawa couplings, we get an 
analytic formula for tan(023)- The phases 4>2,4'3, appearing in Eqs.( |A.9|A.10l ), are fixed in order 



to have real tan(023)- The next step, as shown in Eq. (A. 2), is to apply the rotation i?i3 which 
diagonalises the outer 13 block. Similar to the previous step, this rotation modifies the matrix 
by putting zeros in the 13, 31 entries. 

After applying the 13 rotation, the neutrino mass matrix can be written as, 

mil mi2 

R'ifR2fP2''*mLLP2^R23Ru = | ^12 ^22 | (A.ll) 

m'g 

The last step of the diagonalisation involves modifying the upper 12 block of the matrix. To do 
this, we first multiply the result of the last step by the phase matrix Pi, which introduces the 
phase X- We then apply the real rotation P12. The neutrino mass matrix can then be written 
as follows, 

■ mil rh[^ \ / m; 

P^s^ I m'la m'22 \r'(2=\ m'^ | , (A. 12) 

m'o / \ mi 
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where the new masses 171^2, ^^22 written as, 



mi2 



mi2 e 



rh'22 = m22 e*^^. 



From Eq. (A.12), we can see that the neutrino matrix is successfully diagonalised, however, 



we still need to multiply the result by the phase matrix P3 in order to make all the diagonal 
elements real. To proceed, we write the resulting mass matrix by substituting for the diagonal 
mass terms as m'^ = rme^^i, i = {1, 2, 3}. We then apply the phase matrix and write the phases 
Ui as Ui = (f>i/2. These phases cancel with the phases of the neutrino mass matrix which gives 
a real, diagonal, neutrino matrix as required. 

B Derivation of mass terms 

In this Appendix, we present the derivations of the mass terms resulting from the diagonalisation 



of the mass matrix. After applying the rotation R23 (c.f. Eq. (A.3)), we can derive expressions 



for the masses m'^ and 77122 which are necessary for deriving expressions for ^13 and tan(0i2). To 



find these masses, we first take the trace of both sides of Eq. (A. 10) which gives. 



m22 + rnn 



m22e 
e 



-i2</.3 



Y 



(l + e3(^^3)' + (4)'(e2 + ^l)), 



(B.l) 



We can also express the determinant of both sides of Eq. (|A.10|). This reads, 

m22e" 



^22^3 



2^'^2^33e-2*'^3 - (m23)=^e-2^('^2+'^3) 

/2 52 „2 thf_^^\2 



e", (bf-ec) 



XY 



(B.2) 



We take the mass term m'^ to have the following form, 



(B.3) 



where the parameter f3 is given by Eq. (4.6) and the mass term 771,3 = ''^3lei=o,r?i=o is given by. 



m 



/O ^ j(2</.e-202) 



\e\' + \f\' 
Y 



Using Eqs.( |B.l[ ), ( |B.2[ ), ( |B.4[ ), 77122 can be written as, 

,X ^ Y _ 



77722 « { m^2 + ^2 ( ^ + ^ ) ) (1 - /3) 



where 



™0 _ ~ I ^ ^-2#2 V"'-23 
"1-22 = "7.22|c'=0,ei=0 ~ e ^ 



X 



(B.4) 



(B.5) 



(B.6) 



and the parameter 7/2 is given by Eq. (4.16). 
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In addition to the derivation of the masses 17122, "m'^, applying the rotation R23, modifies the 
masses mi2,TOi3. These become mi2,?7ii3, after diagonahsing the 23 block, and can be derived 
to leading order as follows 



1^12 = C23"ii2e "^^ - S23mi3e 

- -?2(l-7(4)^)-e-^<^3^4(| + f), 



(B.7) 



mi3 



S23"ii2e + C23ini3e 
<3i^-li43f) + e-"'^l4,i"^ + y), 



(B. 



where the parameter 7 is given by Eq. (4.4), the masses m'j^2 = ^12 
mi3|c'=o,ei=o are given to leading order, as presented in [5j, by 



c'=0,ei=0 



and m' 



^0 



13 



m 



12 



O^ggi(</.e-0/)) 



X 



m 



13 



-i<f>2 



X 



S236 + C9c.ce 



,0 ^.J{<l>e 



23'^ 



Y 



(B.9) 
(B.IO) 



After applying the i?i3 rotation, we obtain another mass term, mn, which can be presented 
to leading order as 



m 



mil 



m-ii 



13 



m?i(l - 27(4)' - /3) + ^(27(4)' + /?) + ^(27(4)' + /?) 
2e-^^^^e6(4& + 4ce^^^^-^^^)74 f ^ + 



\X Y 

where the leading order form of mJi = mii|^.=o,ei=o is given in [5], 



2d 



m 



a(46 + c03ce*(<^-^/)) 
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The small parameter eg is written as, 

12 , I f|2N -1 



(B.ll) 



(B.12) 



£6 



X 



\e\' + \ f\' 
Y 



Similar to the derivation of the masses 77122 , ^n'^ , the neutrino masses m'^ and m'2 can be 
written using the trace and the determinant of the upper 12 block of Eq. (A. 12). The real 
neutrino masses mi,m2 can then be written, in the SD cases, as 



mi 



1712 



X 



X|e|2 



+ 



3X' 



1 



Y\bl_\d] 
X\eP Id 



cos(2(/)c') 



The neutrino mass can be written in the SD cases, using Eqs.(B.3), (B.4), as 

'2|e|2 



7713 



Y 



Ic'P \ 
+ ^cos(2(/)c/) j 



(B.13) 
(B.14) 

(B.15) 
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